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Abstract. To establish the electron energy distribution function (EEDF), the second 
derivative of a Langmuir probe current-voltage (I-V) characteristic is numerically 
integrated using the Tikhonov singular value decomposition regularized method. A 
comparison of the numerically intagrated EDDF and by a least-squares fitting is 
discussed. The used I-V characteristic is measured in an ECR plasma source using a 
cylindrical probe and the plasma parameters are determined by the Laframboisc theory. 
This technique allows a rapid analysis of plasma parameters at any gas pressure. The 
obtained EEDF, for the case of the ECR plasma source, shows the existence of two 
groups of electrons with different temperatures. This result is associated with the 
collisional mechanism heating taking place in ECR plasma sources, where low pressure 
plasma is sustained by electron impact ionization of the ground state molecules or 
atoms by energetic electrons arising in the resonance zone. 
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1. Introduction. 

The Langmuir probe is one of the simplest diagnostics tools of the study of ionized 
gases. The probe usually consists of a small sphere or a circular cylinder which is placed 
in the plasma at the point of interest. Some external circuitry is provided so that the 
electric potential of the probe can be varied. A plot of the total current flowing between 
the probe and the plasma versus the probe potential is called the current-voltage (I- 
V) characteristic. Langmuir probes are highly effective in determining the electron 
energy distribution function (EEDF), allowing both spatial and temporal resolution. 
The EEDF is given by the second derivative of the probe I-V characteristic. 

In recent years, several groups have successfully used the numerical differentiation of 
digital probe trace data. [H 121 El- Yet, numerically, the second derivative is a hard task 
and is classified as an ill-posed problem. Ill-posed problems do not have the properties 
of existence, uniqueness and stability e 4j. One of the alternatives to solve this problem 
is the introduction of the concept of conditional well-posed problems [I] . On the other 
hand, the problem of integration of the I-V characteristic has evolved along with the 
solution of the Fredholm integral equation of first class, which is obtained from the 
integration of a second order differential equation. One of the most stable methods to 
solve the Fredholm integral equation are the regularization methods. 

In order to characterize the processes occurring in plasmas it is advisable to obtain 
the electron energy distribution function with a high degree of accuracy. In particular, 
there has been much interest in the electron cyclotron resonance (ECR) plasma sources 
motivated by their applications in ionized physical vapor deposition techniques (I-PVD) 

In this paper, the second derivative of a Langmuir probe current-voltage (I-V) 
characteristic is numerically integrated using the Tikhonov singular value decomposition 
(SVD) regularized method in order to establish the EEDF as the most stable method. 
The existence of two groups of electrons with different temperatures, where low pressure 
plasma is sustained by electron impact ionization of the ground state molecules or atoms 
by energetic electrons arising in the resonance zone is shown. Also, it is important 
to notice that the rate of convergence of the Tikhonov method is much faster than 
that achieved with the least-squares with an orthogonal decomposition (LSQR) and 
truncated singular value decomposition (TSDV) methods. A comparison between the 
EEDF obtained by the Tikhonov regularization method and the one resulting from a 
least-squares fitting is discussed. 

2. Basic equations. 

The electric potential profile near the probe can influence its current collection by setting 
up "barriers of effective potential" around the probe, thereby preventing some particles 
from reaching it. This profile can be strongly influenced by the space charge of the 
particles themselves. However, when r p /\£> is small enough (the sheath around the probe 
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is comparatively thick), this barriers disappears, and the current collection is maximized. 
The currents collected under these conditions then become the orbit-limited-currents. 
This assumption is backed by the orbital motion limit (OML) theory |Hj which implies a 
thick, collisionless sheath (A e r p ). Thus, the calculation of Laframboise P, based on 
the more complete theory of Berstein and Rabinowitz [Jj, showed that the OML limit 
is achived for A^ ~ r p . Here r p is the probe radius (spherical or cylindrical), A^ is the 
Debye radius, and A e is the electron mean free path, respectively. 

When we assume the presence of a stopping field (A e > r p > Ad), as demonstrated 
by Druyvesteyn in [HI EI, the second derivative of the current I with respect to the 
voltage V of the probe is proportional to the electron energy distribution function as 
follows 



where the energy is e = e(V p — V s ), V p is the probe voltage, V s is the plasma potential, 
I is the electron current to the probe, S is the probe surface area, and m, e are the 
electron mass and charge, respectively. The contribution of the ion current to the 
second derivative of the probe current is normally insignificant, and it is therefore not 
subtracted in these EEDF calculations. 
The equation (jlj is in the form 

al 2 I 



dV* 



F(V P ), (2) 



where F(V P ) is a function of the potential on the probe V p with the boundary conditions 
for the current given by 

I\ Vp =a = /(a), 

I\v p = b =I{b). (3) 
A double integration in V p of equation (J2J) results in jTUj 

h(v p ) = f k (v p , v;) f (v;) dv;, (4) 



where 



h (y p ) = i (v p ) - h iy p ) , 

l b-a ' 11 V P ^ V P 



and, 



h(V p ) = 3- [1(b) - 1(a)] + J- [67(a) - al(b)} . (6) 

This integral equation is known as the Fredholm equation of first class and it can be 
solved numerically by regularization methods. The Galerkin discretization becomes the 
more efficient scheme for this class of problems [TT] IT2] . 
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Table 1. Values of constants appearing in H20JI and (|21|1 obtained by a least-squares 
fitting. 
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3. EEDF by the Tikhonov regularization 
method. 

In order to solve the Fredholm's first class equation (0J), the kernel K can be expressed 
in the form of an expansion in terms of singular functions Ui(x), Vi(x), and singular 
values ki. 

n 

K ( x , V)=J2 hui(x)vi(y). (7) 
i=i 

In the occurrence of symmetry, this reduces itself to an eigenfunction expansion. It is 
not necessary that the expansion should converge pointwise; all we require for an L 2 
kernel is that 

ffin J J i^K(x,y) -^2kiUi(x)Vi(y)^ dxdy = 0, (8) 

for which it is necessary and sufficient that the series J2i kf converge. 

When we have an ill-conditioned matrix K, the ki values usually becomes very 
small. In this case the solution of the system (JJJ) proposed by Tikhonov consists in 
replacing the ill-posed problem with a stable minimization problem involving a small 
positive parameter a: instead of attempting to solve the equation (J1J) directly, we seek 
to minimize the quadratic functional [3] 
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Kf-h +a\\Lf\\\K er>> n, (9) 

where L is some linear operator and h denotes h + Sh. If L is suitably chosen, then the 
second term has a smoothing or stabilizing effect on the solution. We may, for example, 
take Lf = f, /', or if the k — th derivative is selected, the process is termed k — th 
order regularization. In the case Lf = f, the solution of the minimization problem (|10|) 
is then obtained as the solution of the linear equation 

(K T K + al) f a = K T h. (10) 

The operator acting on f a is clearly positive-definite when a > and consequently it 
has a bounded inverse. Solving in terms of singular functions we obtain 

i Kj -+- Ct 

Comparing this with the exact expansion we can see that the effect of regularization 
has been to insert a filter factor ki/{kf + a). This is close to unity so long as ki is large 
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Table 2. Values of plasma parameters obtained by the Langmuir theory, the Tikhonov 
regularization method and by a least-squares fitting. 
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compared with a but tends to zero as k{ — > 0, the rate of transition depending on a. If 
we split h into h + Sh the expression (jllj) becomes 

/^E^^H^^G/)- (12) 

As for the first term, it is advantageous to make a small in order to reduce the error due 
to regularization; by contrast, the second term, which only consists of error, is made 
small by taking a large. Thus, there is a conflict, and we would like to achieve the best 
compromise jllj . 

4. Least-squares fitting. 

In the OML theory of ion collection, the ion current flowing to a negatively biased probe 
is independent of the shape of the plasma potential V(r) as long as the current is limited 
only by the angular momentum of the orbiting ions [Hj. This requires the arbitrary 
assumption of either a "sheath edge" s, beyond which the ion energy distribution is 
Maxwellian, or a V(r) varying so slowly that no "absorption radius", inside of which 
all ions are drawn in, exists between the probe and infinity This condition is never 
satisfied even at modest densities. For s — ► oo and a Maxwellian ion distribution at 
temperature Tj, the OML current to a cylinder probe is given by 

2 



Sjr 



X 1/2 + e* (l - eri( X 1/2 ~ 



Sj r ^=y/T+^, (13) 



where x = —eV p /kTi and j r is the random thermal ion current. As Tj — > 0, the Tj 
dependencies of x an d jr vanish, and a finite limiting value of the OML current is 
reached [TH] 

In the OML theory, the ion current of equation (JT4"jl can be represented as 

i{v p ) = -p(v p -v f y + c , (is) 
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where (3, 7, Vf and C are constants to be fitted. In our case Vf acquires the meaning 
of the float potential. 

The three trial functions in the least-squares fitting of the I-V characteristic for the 
electronic component have the form 

F 1 (V P ) = AiiVp - V 8 )e^^- V '\ (16) 

and 

K(V p ) = A i (V p -V,) 2 e-Wp-v.)]\ ( i = 2)3 ) (17) 

where A4, Ui, (i = 2, 3) and V s are constants that must be valued. Here V s adopts the 
character of the plasma potential. 

Substituting equation (JTfij) into an expression for the current to the probe, written 
in the form [TH] 

poo 

I{V) = B X (V p -x)F(V p )dV p , (18) 

J X 

where B\ = (4/eS')y / l/2m e , x = eV p and in place of F(V P ) we take any of functions 
(|16 |) -([17 |) . we obtain for the electronic component 



h(V p ) A x 



[2 + Ul (V p - V s )] exp l- Ul (V p - V s )] , (19) 



Bi ui 

where the subindex 1 in / refers to trial function ([16)1 . Analogously, from relations (|T7jl 
we get 

Ij(V P ) _ Aj 



^2 

' s 



2 {eexp [-u 2 (V p - V s ] 

i 

+ u i ^{Vp-V s ){^f[u i {Vp-V s )]-l)} (i = 2,3) (20) 



5. Analysis and discussion. 

With a view to obtaining accurate values of n, T e and V s (but not Tj) from the EEDF, 
we shall illustrate a procedure using data obtained in an ECR discharge with P = 10 
mTorr in argon gas, taken with an rf-compensated cylindrical probe with a radius 0.4 
mm and 4 mm in length. 

It is considered that the current is collected within the area at the tip of the probe 
expressed by S = 2ttR p L + ttR 2 , where L is the probe length. The entire I-V curve in 
this example is shown in figure H 

From the Langmuir theory, assuming a Maxwellian distribution for electrons |llHj . 
we obtain in the case of the I — V curve of figure [TJ that T e = 6.88 eV, V s = 22.078 
V, and n e = 1.39 x 10 16 m~ 3 . Thus, the first derivative shown in figure ^ as well as 
experimental data are used. Here, it is important to notice that the extremum of the 
first derivative is not completely well defined, as observed in figure H 

The EEDF is calculated following the procedure described in §3, by the integration 
of equation (0J) with kernel (0), and using the Galerkin discretization method |14j . 
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Table 3. Values of plasma parameters obtained by the Langmuir theory, the Tikhonov 
regularization method and by a least-squares fitting for the energetic group of electrons. 
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In figure |3]the EEDF obtained by the Tikhonov technique of zeroth order is shown, 
along with the corresponding values of the regularization parameters a = (0.09, 0.08, 
0.07, 0.06, 0.05, 0.04, 0.03, 0.02). 

From these charts we gather that the most probable energy is T e = 5.66 eV. With 
this electron temperature, it is easy to calculate the plasma electron density by the 
relation 

1 /2 

'ft'^-^y ■ (21) 

obtaining that n e = 1.54 x 10 16 m -3 . It is also important to notice that the 
maximum value of the EEDF by the Tikhonov regularization method takes the value 
of F(e)\ max = 7.34 x 10 14 m~ 3 eV -1 . After an analysis of several calculations, we can 
say that the bumps observed in the EEDF are related with numerical effects and are 
originated by the irregularities of the I-V curve. 

These results are validated in the following by a least squares fitting. As it has 
been described in §IV about the ionic part of the I-V characteristic, we fit the data to 
equation (fTB^) determining the values for the constants 7, (3 and C . These values are 
summarized in Table ^ This fitting is shown in figure ^ In figure El the measured 
values I 2 versus V curves are plotted ^3]. From these charts, we can observe that the 
obtained values for the constants show a good agreement with the data values. 

Analogously, for the electronic component in the I-V curve, we now employ 
expressions (fL8|) and (p"9|) for a least squares fit to the data. The corresponding 
values of constants appearing in (fTH|) and (fT9*|) are given in Table [D Then, we can 
express the EEDF as a sum of functions (|16 p - (fTTj) . This distribution is plotted in 
figure E Introducing the corresponding dimensions, we obtain a comparison with the 
EEDF obtained by the Tikhonov regularization method shown in figure El For the 
energetic group of electrons we obtained that T e = 18.831 eV, n e = 8.66-E15 m~ 3 , and 
F(e)\ max = 3.42.ET3 m -3 eV -1 . From the least-squares distribution we obtain for the 
main group of electrons that the most probable energy is T e = 5.331 eV. The maximum 
value found for the EEDF is F(e)\ max = 1.131 x 10 15 m~ 3 eV -1 . Considering this value 
for the temperature T e , we find by equation (|2~Tj) the density n e = 1.62 x 10 16 m~ 3 . 
For the values computed after fitting, it is possible to accomplish the integration of the 
adjusted characteristic, to obtain [13 HE] 

POO 

n e = / F(e)de = 7.8 x 10 15 m" 3 , (22) 
Jo 
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and 

T eff = -(e) = 2(3<T 1 / sF(e)de = 8.40 eV. (23) 
o Jo 

Here, the related problem of the indetermination of the zero position in the EEDF, 
as it is observed from figure H3 deserves some attention. This problem will always be 
present when no electronic saturation is reached. Several results obtained by applying 
the two methods are summarized in Tableland Table El 
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Figure 1. Sample I-V curve to be analyzed obtained from an 0.4 mm diam, 4 mm 
long probe in an ECR discharge in 10 mTorr argon gas (solid). First derivative of the 
FV curve (dotted), and saturation ion and electronic currents computed after fitting 
(dashed). 



Figure 2. Square of saturation ion current vs probe voltage as measured (solid) and 
as computed after fitting (dashed). 



Figure 3. EEDF resulting from the least squares fitting (solid) and by the Tikhonov 
regularization method (dotted). Here are shone the locations of the main and energetic 
groups of electrons respecting the energy. 
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